Method for Estimating Relative Clock Frequency Offsets to Improve Radio Ranging Errors

ABSTRACT

A method and system improves two-way radio ranging accuracy by estimating a relative clock frequency offset between a first clock X of a first transceiver and a second clock Y a second transceiver. The first transceiver transmits a first packet at time t 0  received by the second transceiver at a time t 1 . The second transceiver transmits a second packet at a time t 2  received by first transceiver at a time t 3 . The second transceiver transmits a third packet at a time t 4  received at a time t 5 . The relative clock frequency offset is then 
     
       
         
           
             
               
                 Δ 
                  
                 
                     
                 
                  
                 
                   f 
                   XY 
                   ′ 
                 
               
               ≅ 
               
                 
                   2 
                    
                   
                       
                   
                    
                   
                     f 
                      
                     
                       ( 
                       
                         
                           N 
                           24 
                           Y 
                         
                         - 
                         
                           N 
                           35 
                           X 
                         
                       
                       ) 
                     
                   
                 
                 
                   
                     N 
                     24 
                     Y 
                   
                   + 
                   
                     N 
                     35 
                     X 
                   
                 
               
             
             , 
           
         
       
     
     where f is a nominal clock frequency of the first and second clocks, N Y   24  is a measured first delay between times t 2  and t 4  of the second clock, N X   35  is a measured second delay between times t 3  and t 5  of the first clock X.

FIELD OF THE INVENTION

This invention relates to wireless communication systems, and more particularly, to accurate ranging estimation through relative clock frequency offsets (RCM) compensation between wireless devices.

BACKGROUND OF THE INVENTION

In many wireless communication networks, devices are not synchronized to a single reference clock. Instead, each device relies on its own clock. Due to the limitation of the hardware, manufacturing processes, temperature drift, component aging, etc., the actual operating frequency of the clock is generally different from the designed nominal frequency. For example, if a 100 MHz clock signal is generated by a crystal that has a tolerance of +/−20 parts per million (ppm), the actual frequency of the clock can be any value in a range from 99,998,000 Hz to 100,002,000 Hz.

A difference between the actual frequency of the clock and the designed nominal frequency is referred to as the ‘clock frequency offset’ (CFO) or ‘absolute clock frequency offset’ (ACFO). The difference between the frequencies of two independent clocks (with the same nominal frequency) is called ‘relative clock frequency offset’ (RCFO). For example, for two clocks with nominal frequency of 100 MHz, if one clock has a frequency of 99,998,000 Hz and the other clock has a frequency of 100,002,000 Hz. The RCFO is 4 KHz, or 40 ppm.

The CFO and RCFO can cause problem. When the clock is used for time measurement, the ACFO and/or RCFO introduce errors. In digital circuit, time is generally measured as the number of clock cycles between two time instances

${N = {\frac{t}{T_{actual}} = {t \cdot f_{actual}}}},$

where T_(actual) is the actual clock period, and f_(actual) is the actual clock frequency. As an example, a timer based on an accurate 100 MHz clock will count 100,000 clock cycles in a 1 millisecond period. If the ACM is −20 ppm, the actual count is ˜99,998 clock cycles instead. If the ACFO is not compensated, which is often the case when the device has not knowledge of the ACFO of its own clock, then the time is simply computed using the nominal frequency:

$\begin{matrix} \begin{matrix} {t = {N \cdot T_{norm}}} \\ {= \frac{N}{f_{norm}}} \\ {= {t \cdot \frac{f_{actual}}{f_{norm}}}} \\ {= {t + {t \cdot \frac{\Delta \; f_{abs}}{f_{norm}}}}} \end{matrix} & (1) \end{matrix}$

where T_(norm) is the nominal clock period, f_(norm) is the nominal clock frequency, and Δf_(abs) is the ACFO. In the example, f_(norm)=100 MHz, t=999,980 ns, which results in a 20 ns time error.

If the same time period is measured by two devices using independent clocks, the RCFO causes a discrepancy in the measurements. For example, if the nominal clock frequency is 100 MHz for both timer X and timer Y, timer X has a −20 ppm offset and timer Y has a +20 ppm offset, which is a RCFO of 40 ppm. A 1 msec time period will be measured as 99,998,000 ns and 100,002,000 ns respectively for timer X and timer Y, a difference of 40 ns.

The RCFO has a significant impact on the ranging accuracy of a two-way TOA (TW-TOA) system. A TW-TOA is a method can also be used for ranging between two transceivers. TW-TOA does not require exact synchronization between the transmitter and receiver. The TW-TOA ranging method is used in the IEEE 802.15.4a standard.

In TW-TOA, two transceivers exchange packets and the round trip delay is measured. A typical exchange is as follows: the first device transmits a packet to the second device. After receiving the packet, the second device transmits a packet back to the first device. The first device measures the total time elapsed from the transmission of the first packet to the reception of the second packet. The second device measures the time elapsed from the reception of the first packet to the transmission of the second packet. The round trip travel time of the signal is calculated as the difference between these two measured values. The distance between the devices is calculated as the product of one half of the round trip traveling time multiplied by the speed of the signal (3×10⁸ meter/sec in free space for electromagnetic signal).

Because the turnaround time at the second device is much larger than the signal ‘flying’ time, the error caused by RCFO is the dominant factor in the overall time estimation error for a TW-TOA system. If the devices can obtain an accurate estimation of the RCFO, then the time estimation accuracy can be improved by compensating for it. Clearly, an accurate RCFO estimation is very important in achieving accurate time measurement.

FIG. 4 shows a prior art ranging procedure. Device X transmits a request packet (ranging-req) 111 at time t₀. The packet is received by device Y at time t₁. After some delay, device Y transmits an acknowledge packet (ranging-ack) 112 at time t₂. The range-ack packet is received by device X at time t₃. The turn-around delay t₁₂ is measured by device Y and sent to device X either in the ranging-ack packet, or in a separate packet (ranging-data). The delay t₀₃ is measured by device X.

Let t_(p) be the ‘fly’ time of the signal. Without considering the frequency offsets, the estimated round-trip fly time is

2.t′ _(p) =t ^(X) ₀₃ −t ^(Y) ₁₂,

and t^(X) ₀₃ the measured round-trip time from the transmission of the range-request packet to the reception of the range-acknowledge packet by device X based on its clock clk_(X), and t^(Y) ₁₂ is the measured turn-around time from the beginning of the received ranging-request packet to the beginning of the transmitted ranging-acknowledge packet by device Y based on its clock clk_(Y).

The range between device X and Y is

D_(XY)=t′_(p).C,

where C is the speed of the electromagnetic signals, e.g., 3×10⁸ in/sec in free space.

Because the true round-trip fly time is

2.t _(p) =t ₀₃ −t ₁₂,

t ^(X) ₀₃ =N ^(X) ₀₃ /f=t ₀₃ .f _(X) /f, and

t ^(Y) ₁₂ =n ^(Y) ₁₂ /f=t ₁₂ .f _(Y) /f, and

the estimation error is

Δt _(p) =t′ _(p) −t _(p)=(t ^(X) ₀₃ −t ^(Y) ₁₂)−(t ₀₃ −t ₁₂).

Given that the process time is much longer than the fly time and t₁₂>>2t_(p),

Δt _(p)=(t ^(X) ₀₃ −t ₀₃)−(t ^(Y) ₁₂ −t ₁₂)=t ₁₂(Δf _(XY) /f)+2t _(p)(Δf _(X) /f≈t ₁₂(Δf _(XY) /f),

where Δf_(X) is the frequency offset of the timer clock with respect to the nominal frequency f and Δf_(XY) is the relative frequency offset between the timer clock for device X and the timer clock for device Y. The above equation indicates that reducing t₁₂ and Δf_(XY)/f will improve the time estimation accuracy. However, in practice, the desired improvements may not be feasible. Δf_(XY)/f is determined by the components performance (in particular, the crystal oscillator) and a higher accuracy/stability components are more costly and may require more complicated circuit design.

The process time, t₁₂ is often determined by several factors and cannot be easily reduced because many communication standards impose a quite period between transmissions, and the receiving device must process the data and pass the packet to the upper layer and the processing time is long, and the length of the payload portion of the packet is often non-zero.

Conventional methods for estimating the RCF(in wireless systems include the following.

One method uses the preamble of a received packet. Generally the preamble of the packet includes multiple symbols. By measuring a length of the preamble or a time interval between two symbols, the receiving device can estimate the RCFO of its clock with respect to the clock of the transmitting device. Such a method does not yield an accurate RCFO estimation because the preamble has a finite duration. Therefore, the time interval between two symbols cannot be very large, the method does not benefit from the processing gain by using as many symbols in the preamble, and the method is sensitive to carrier frequency offset and sample timing errors.

A symmetric double-sided two-way protocol exchange data/ranging packet between two devices in the following sequence:

-   -   Device A transmits a first packet;     -   Device B receives the first packet;     -   Device B transmits a second packet;     -   Device A receives the second packet;     -   Device A transmits a third packet; and     -   Device B receives the third packet.     -   Devices A and B exchange the timing information in additional         data packets.

Device A measures a round trip time t_(roundA) between transmitting the first packet and receiving the second packet. Device B measures a reply time t_(replyB) between receiving the first packet and transmitting the second packet, and device A measures a reply time t_(replyA) between receiving the second packet and transmitting the third packet.

The drawbacks of such a method are: the times t_(replyA) and t_(replyB) must be identical or very close in order to cancel out the error caused by the RCFO; additional transmissions are needed to pass the information from one device to the other device; and delay mismatches in RF frond ends of two devices may introduce additional errors.

Another method uses a combination of preamble and payload section or possibly post amble, see U.S. patent application Ser. No. 11/749,517, “Method for Reducing Radio Ranging Errors Due to Clock Frequency Offsets,” filed by Sahinoglu on Jun. 30, 2007, incorporated here in by reference.

Another method uses frequency domain analysis. By performing a fast Fourier transform (FFT) on part or all of the preamble of a packet, the RCFO can be estimated. This approach requires the storage of a large number of samples, a complex FFT operation and therefore requires additional hardware, software and power. The accuracy is dependent on the signal-to-noise ratio (SNR) of the received signal and the size of the preamble portion used for FFT.

It is desired to provide a method for estimating the RCFO without the above problems and complexities.

SUMMARY OF THE INVENTION

The invention provides a method for estimating a relative clock frequency offset (RCFQ) between two wireless communication devices with asynchronous clocks by performing a data packet exchange. The method minimizes range estimation error in a two-way time-of-arrival (TOA) measurements.

A first transceiver transmits a first packet. A second transceiver transmits a second packet, and then, after a first delay a third packet. The first transceiver receives the second packet, and then after second delay the third packet. The relative clock frequency offset is determined from the measured first and second delays.

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1 is a block diagram of two devices communicating with each other by exchanging packets according to an embodiment of the invention.

FIG. 2 is a block diagram of a timer according to an embodiment of the invention;

FIG. 3 is a timing diagram of time error caused by clock frequency offsets; and.

FIG. 4 is a timing diagram of a packet exchanges protocol between two devices without RCFO estimation; and

FIG. 5 is a timing diagram of a packet exchanges protocol between two devices according to an embodiment of the invention.

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS

Wireless Ranging Devices

FIG. 1 shows wireless communication device X 101 and device Y 102 according to an embodiment of the invention. A distance between device X and device Y is D_(XY) 110. Device X can transmit a data packet 120 to device Y, and device Y can transmit a data packet 130 to device Y. FIG. 1 also shows timer X 141 and timer Y 142, and clock signals clk_(X) 151 and clk_(Y) 152. The timers count the clock signals. The timers are identical other than that count clock signals having different frequencies due to offsets from their nominal (design) frequency. A device may have multiple timers.

Timers

FIG. 2 generally shows a timer 201 operated according to a signal clk 202 having a nominal clock frequency f_(norm). Of course, the nominal frequency is known at both devices. The time 201 forms the basis for timers 141-142.

The timer starts counting when a signal cnt_(start) 211 is active (ON), and stops counting when cnt_(stop) 212 is active. The timer outputs a count N 213 and the time elapsed from cnt_(start) to cnt_(stop) can be determined using

N=(t _(stop) −t _(start))f _(norm).

Relative Clock Frequency Offsets (RFCO)

FIG. 3 shows a time measurement error caused by clock frequency offset. At a time t₀, the cnt_(start) signal is active, and the timer starts counting. Then, at a time t₁, the cnt_(stop) is active and the timer stops counting.

If the nominal clock frequency of the timer is f_(norm), then the output count is N. If the actual clock frequency of timer X is f_(X), then the output count of timer X is N_(X). If the actual clock frequency of time Y is f_(Y), then the output count of time Y is N_(Y) given as follows:

N=t·f;

N _(X) =t·f _(X) =t·(f+Δf _(X))=N+t·Δf _(X) =N+ΔN _(X);

N _(Y) =t·f _(Y) =t·(f+Δf _(Y))=N+t·Δf _(Y) =N+ΔN _(Y); and

ΔN=N _(X) −N _(Y)=(N+t·Δf _(X))−(N+t·Δf _(Y))=t·(Δf _(X) −Δf _(Y))=t·Δf _(XY),

where Δf_(X) is the frequency error of clock X with respect to the nominal frequency f; Δf_(Y) is the frequency error of clock Y with respect to the nominal frequency f; Δf_(XY) is the relative frequency error between clock X and clock Y. ΔN_(X) and ΔN_(Y) are the count error with, respect to the nominal count N.

If N_(avg)=(N_(X)+N_(Y))/2, then t can be approximated as:

${t = {{\frac{\left( {N_{X} + N_{Y}} \right)}{2\; f} - \frac{\left( {{\Delta \; N_{X}} + {\Delta \; N_{Y}}} \right)}{2\; f}} \approx \frac{N_{avg}}{f}}},{and}$ ${\Delta \; N} = {{{{t \cdot \Delta}\; f_{XY}} \approx {{\frac{N_{avg}}{f} \cdot \Delta}\; f_{XY}}} = {{\frac{N_{X} + N_{Y}}{2 \cdot f} \cdot \Delta}\; {f_{XY}.}}}$

If N_(X), N_(Y) and the nominal frequency f are known, the relative clock frequency offset is

${\Delta \; f_{XY}} \cong {\frac{2\; {f\left( {N_{X} - N_{Y}} \right)}}{N_{X} + N_{Y}}.}$

RCFO Protocol

FIG. 5 shows a packet exchange protocol and an associated RCF( estimation method according to an embodiment of the invention. In FIG. 5, the time axis runs down. In the following, the timing information at device X is with respect to clock X and timer 151, and at device Y with respect to clock Y and time 152. This is implicitly understood and explicitly indicated in the superscripts and subscripts. Thus, in FIG. 5, the times t are equivalent to the clock counts N, as shown.

At time t₀, device X transmits a ranging request packet 111. Device Y receives the ranging request packet at time t₁.

After receiving the ranging request packet and a delay of t₁₂, device Y transmits a ranging acknowledgement packet 112 at time t₂. Device X receives the ranging acknowledgement packet at time t₃.

After receiving the ranging request packet. 111 and a first delay t₂₄, device Y transmits a ranging data 113 at time t₄. Device X receives the ranging data packet at time t₅.

The estimated round trip time without considering the frequency offsets is

2·{dot over (t)} _(p) =t ^(X) ₀₃ −t ^(Y) ₁₂,

where, t^(X) ₀₃ is the measured round trip time transmitting the ranging request packet to receiving the range acknowledge packet by device X determined as

${t_{03}^{X} = {\frac{N_{03}^{X}}{f} = \frac{t_{03} \cdot f_{X}}{f}}},$

and t^(Y) ₁₂ is the measured turn-around time from receiving the ranging request packet to transmitting the ranging acknowledge packet by device Y determined as

$t_{12}^{Y} = {\frac{N_{12}^{Y}}{f} = {\frac{t_{12} \cdot f_{Y}}{f}.}}$

The range or distance 110 between device X and Y is

D _(XY) =t′ _(p) ·C

where C is the speed of light.

Because the true round trip time is

2·t _(p) =t ₀₃ −t ₁₂,

the estimation error is

Δt _(p) ={dot over (t)} _(p) −t _(p)=(t ^(X) ₀₃ −t ^(Y) ₁₂)−(t ₀₃ −t ₁₂).

Given that the processing time is much longer than the signal ‘flying’ time, i.e.,

$\begin{matrix} {{{t_{12}\operatorname{>>}{2\; t_{p}}},{then}}\begin{matrix} {{\Delta \; t_{p}} = {\left( {t_{03}^{X} - t_{03}} \right) - \left( {t_{12}^{Y} - t_{12}} \right)}} \\ {= {{t_{12}\left( \frac{\Delta \; f_{XY}}{f} \right)} + {2{t_{p}\left( \frac{\Delta \; f_{X}}{f} \right)}}}} \\ {{\approx {t_{12}\left( \frac{\Delta \; f_{XY}}{f} \right)}},} \end{matrix}} & (1) \end{matrix}$

where Δf_(X) is the frequency offset of the timer clock with respect to the nominal frequency f and Δf_(XY) is the relative frequency offset between the timer clock for device X and the timer clock for device Y. Equation (1) indicates that reducing t₁₂ and Δf_(XY)/f improves the time estimation accuracy.

However, in practice, the desired improvements may not be feasible. The value Δf_(XY)/f is determined by the components performance, in particular, the crystal oscillator, and a higher accuracy/stability components are more costly and may require more complicated circuitry design.

The processing time t₁₂ is often determined by several factors and cannot be reduced easily. These factors include the following. Communication standards usually impose a ‘quite’ period between transmissions. The receiving device must process the data and pass the packet to an upper layer and the processing time is substantially longer than the round trip time. The packets have a non-zero length.

However, if we obtain an accurate estimation of Δf_(XY), then we can improve the accuracy of the time estimation. The flight time estimation of the signals of each packet between the devices (transceivers), with relative clock frequency offset, is

$\begin{matrix} {{2 \cdot t_{p}^{C}} = {t_{03}^{X} - {t_{12}^{Y} \cdot {\left( {1 - \frac{\Delta \; f_{XY}^{\prime}}{f}} \right).}}}} & (2) \end{matrix}$

The estimation error is

$\begin{matrix} \begin{matrix} {{\Delta \; t_{p}^{C}} = {{t_{12}\left( \frac{\Delta \; f_{XY}}{f} \right)} - {t_{12}^{Y}\left( \frac{\Delta \; f_{XY}^{\prime}}{f} \right)}}} \\ {= {t_{12}\left\lbrack {\left( \frac{{\Delta \; f_{XY}} - {\Delta \; f_{XY}^{\prime}}}{f} \right) + \left( \frac{\Delta \; {f_{XY} \cdot \Delta}\; f_{XY}^{\prime}}{f^{2}} \right)} \right\rbrack}} \\ {{= {{t_{12}\left( \frac{{\Delta \; f_{XY}} - {\Delta \; f_{XY}^{\prime}}}{f} \right)} + \left( \frac{\Delta \; {f_{XY} \cdot \Delta}\; f_{XY}^{\prime}}{f^{2}} \right)}},} \end{matrix} & (3) \end{matrix}$

where Δf′_(XY) is the estimated RCFO.

Equation (3) shows that an accurate estimation of relative frequency offset reduces the time estimation error significantly. If an accurate relative frequency clock offset estimation is achieved, i.e., Δf′_(XY)=Δf_(XY), then first term in Equation (3) becomes zero. The second term is very small because f>>Δf_(XY) and f>>Δf′_(XY). Therefore, Δt^(C) _(p) is significantly smaller than Δp. Hence, the time estimation is very accurate.

To estimate the RCFO, device X measures the second delay t₃₅, which is the time elapsed between the reception of the ranging acknowledgement packet 112 and the reception of the ranging data packet 113, as counts N^(X) ₃₅. Device Y measures the first delay t₂₄, which is the time elapsed between the reception of the ranging acknowledgement packet 112 and the reception of the ranging data packet 113, as counts N^(Y) ₂₄.

Device Y transmits the time measurement N^(Y) ₂₄ to device X.

Device X determines the relative clock frequency offset (RICO) Δf_(XY) as

$\begin{matrix} {{{\Delta \; f_{XY}^{\prime}} \cong \frac{2\; {f\left( {N_{24}^{Y} - N_{35}^{X}} \right)}}{N_{24}^{Y} + N_{35}^{X}}},} & (3) \end{matrix}$

where f is the nominal frequencies of the first and second clocks.

The first delay N^(Y) ₂₄ is an exact value as it is the delay between two transmitted packets in device Y. The second delay N^(X) ₃₅ is the time difference between two received packets by device X.

Then, device X can then determine t^(C) _(p) using Equation (2).

The value of the first delay N^(Y) ₂₄ can be made known to device X using methods including, but not limited to the following:

-   -   N^(Y) ₂₄ is predetermined by the network;     -   N^(Y) ₂₄ is predetermined by device Y and transmitted in the         ranging acknowledgement packet 112;     -   N^(Y) ₂₄ is transmitted in the ranging data packet 113; and     -   N^(Y) ₂₄ is transmitted to device X in some other packet.

It should be understood, the device Y can also determine the RCFO similarly provided with the second delay.

In the cases when the value of RCFO Δf_(XY) is relatively stable, the RCFO estimation can be performed periodically.

Effect of the Invention

Compared with conventional method, the invention has the following advantages.

The invention does not require a modification of a structure of packets that are used to perform the ranging or data communication.

Any effect of channel impulse response, RF circuit and baseband circuit imperfection are cancelled.

Any quantization error caused by a finite period of the clock can be reduced by increasing the first delay.

The accuracy of estimating the RCM is high.

The requirement for additional hardware, computation and power is very low.

The method is compliant with the IEEE 802.15.4a standard, and does not require a modification to the media access (MAC) protocol used by that standard.

Although the invention has been described by way of examples of preferred embodiments, it is to be understood that various other adaptations and modifications can be made within the spirit and scope of the invention. Therefore, it is the object of the appended claims to cover all such variations and modifications as come within the true spirit and scope of the invention. 

1. A method for improving two-way radio ranging accuracy by estimating a relative clock frequency offset between a first clock X of a first transceiver and a second clock Y a second transceiver, comprising: transmitting, by the first transceiver, a first packet at time t₀ of the first clock X; receiving, by the second transceiver, the first packet at a time t₁ of the second clock Y; transmitting, by the second transceiver, a second packet at a time t₂ of the second clock; receiving, by the first transceiver, the second packet at a time t₃ of the first clock X; transmitting, by the second transceiver, a third packet at a time t₄ of the second clock Y; and receiving, by the first transceiver, the third packet at a time t₅ of the first clock X, and in which the relative clock frequency offset is ${{\Delta \; f_{XY}^{\prime}} \cong \frac{2\; {f\left( {N_{24}^{Y} - N_{35}^{X}} \right)}}{N_{24}^{Y} + N_{35}^{X}}},$ where f is a nominal clock frequency of the first and second clocks, N^(Y) ₂₄ is a measured first delay between times t₂ and t₄ of the second clock, N^(X) ₃₅ is a measured second delay between times t₃ and t₅ of the first clock X.
 2. The method of claim :1, further comprising: measuring a distance between the first transceiver and the second transceiver based on the relative clock frequency offset.
 3. The method of claim 1, in which the first delay is predetermined by a network in which the first and second transceivers operate.
 4. The method of claim 1, in which the first delay is predetermined by the second transceiver and transmitted in the first packet.
 5. The method of claim 1, in which the first delay is transmitted in the second packet.
 6. The method of claim 1, in which the first delay is transmitted by the first transceiver to the second transceiver in a third packet.
 7. The method of claim 1, in which a structure of the first and second packets complies with a communications standard.
 8. The method of claim 1, further comprising: canceling an effect of channel impulse response, RF circuit and baseband circuit imperfection while determining the relative clock frequency offset.
 9. The method of claim 1, further comprising: increasing the first delay to reduce quantization errors.
 10. The method of claim 1, in which an estimated flight-time of the signals between the first and second transceivers is ${t_{p}^{C} = {\frac{1}{2}\left( {t_{03}^{X} - {t_{12}^{Y} \cdot \left( {1 - \frac{\Delta \; f_{XY}^{\prime}}{f}} \right)}} \right)}},$ in which t^(X) ₀₃=t₃−t_(o) measured by the first clock X, and t^(Y) ₁₂=t₂−t₁ measured by clock Y.
 11. A system for improving two-way radio ranging accuracy by estimating a relative clock frequency offset between a first clock X of a first transceiver and a second clock Y a second transceiver, comprising: a first transceiver transmits a first packet at time t₀ of a first clock X; a second transceiver receives the first packet at a time t₁ of a second clock Y, and in which the second transceiver transmits a second packet at a time t₂ of the second clock, which is received by the first transceiver at a time t₃ of the first clock X, and in which the second transceiver transmit a third packet at a time t₄ of the second clock Y, which is received by the first transceiver at a time t₅ of the first clock X, and in which the relative clock frequency offset is ${{\Delta \; f_{XY}^{\prime}} \cong \frac{2\; {f\left( {N_{24}^{Y} - N_{35}^{X}} \right)}}{N_{24}^{Y} + N_{35}^{X}}},$ where f is a nominal clock frequency of the first and second clocks, N^(Y) ₂₄ is a measured first delay between times t₂ and t₄ of the second clock, N^(X) ₃₅ is a measured second delay between times t₃ and t₅ of the first clock X. 